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ABSTRACT 

The D/H ratio in cometary water is believed to be an important indicator of the conditions under which icy planetesimals formed and can provide 
clues to the contribution of comets to the delivery of water and other volatiles to Earth. Available measurements suggest that there is isotopic 
diversity in the comet population. The Herschel Space Observatory revealed an ocean-like ratio in the Jupiter-family comet 103P/Hartley 2, 
whereas most values measured in Oort-cloud comets are twice as high as the ocean D/H ratio. We present here a new measurement of the D/H 
ratio in the water of an Oort-cloud comet. HDO, H 2 0, and H' s O lines were observed with high signal-to-noise ratio in comet C/2009 PI (Garradd) 
using the Herschel HIFI instrument. Spectral maps of two water lines were obtained to constrain the water excitation. The D/H ratio derived from 
the measured H 2 ie O and HDO production rates is (2.06 ± 0.22) x 10~ 4 . This result shows that the D/H in the water of Oort-cloud comets is not as 
high as previously thought, at least for a fraction of the population, hence the paradigm of a single, archetypal D/H ratio for all Oort-cloud comets 
is no longer tenable. Nevertheless, the value measured in C/2009 PI (Garradd) is significantly higher than the Earth's ocean value of 1.558 x 10~ 4 . 
The measured 16 0/ 18 ratio of 523 ± 32 is, however, consistent with the terrestrial value. 

Key words. Comets: general; Comets: individual: C/2009P1 (Garradd); Submillimeter: planetary systems; Oort cloud; Astrochemistry 



1. Introduction 

Having retained and preserved pristine material from the solar 
nebula, comets cont ain unique clues to the history and evolution 
of the solar system dlrvine et all [2000). Isotopic ratios are im- 
portant indicators of the conditions under which cometary ma- 
terials formed, since isotopic fractionation is very sensitive to 
physical conditions. In addition, the characterization of the iso- 
topic composition of long-period comets from the Oort cloud 
and of short-period comets from the Jupiter family can provide 
clues to their formation regions in the early solar system. 

The D/H ratio in water has been determined in several Oort- 
cloud comets (OCC) using different technique s, with most mea- 
surements agreeing with a value of ~ 3 x 10" 4 dJehin et all [2009. 
and references therein). In addition, the ratio wa s measured us- 
ing the ESA Herschel Space Observatory dPilbratt et all [2 010) 
in the Jupiter-family comet (JFC) 103P/Hartley 2 dHartogh et al. . 
1201 lb , where it was found to be (1.61 + 0.24) x 10~ 4 , i.e., con- 
sistent with the Vienna standard mean ocean water (VSMOW) 
value equal to 1.558 x 10~ 4 . This result was unexpected. 
Models quantifying the deuterium enrichment factor in the so- 
lar nebula with respect to the protosolar value predicted an in- 
creas e in the D/H ratio with increasing helioc entric distance 
(e.g.. iHersant et all 120011: IKavelaars et all 1201 lb . According to 
the most accepted theory, OCCs and JFCs originate from the 
same population of objects formed in the Uranus-Neptune zone 
dDones et all 12004). though part of the Oort cloud was possibly 
formed by planetesimals scattered from the Jupiter-Saturn region 



* Herschel is an ESA space observatory with science instruments 
provided by European-led principal investigator consortia and with im- 
portant contribution from NASA. 



(Brasser, 2008). Therefore, JFCs were expected to exhibi t a D/H 
ratio similar to or higher than that in OCCs dKavelaars et all 
1201 ih . 

We present in this paper a new measurement of the D/H 
ratio in the water of an Oort-cloud comet, obtained with the 
same instrumentation and methodology as those used for comet 
103P/Hartley 2. Comet C/2009 PI (Garradd) was observed 
with the Heterodyne Instrument for the Far-Infrared (HIFI, 
Ide Graauw et all l2010h in the framework of the Herschel guar- 
anteed time ke y programme "Water and related chemistry in the 
solar system" (Hart ogh et alll2009l) . 

2. Observations 

Comet C/2009 PI (Garradd) is a long-period comet origi- 
nating from the Oort cloud (P = 127 000 yr ! _orbit incli- 
nation of 106° with respect to the ecliptic; Nakanol 120121) . 
Discovere d on 13 August 2009 at a h eliocentric distance r/ t = 
8.7 AU (McN augh t & Garraddl 12009). it passed perihelion on 
23 December 201 1 at r^ = 1.55 AU. The HDO observations with 
the HIFI instrument were performed on 6 October 2011, when 
the comet was at r/, = 1.88 AU, and a distance from Herschel 
of 1.76 AU. The Hj 8 and H2O lines were observed simultane- 
ously. Since the H?Q ground state rotational lines in comets ar e 
optically thick ( Be nsch & Berginl 120041; IZakharov et all 120071) . 
lines of the rare oxygen isotopic counterpart H^O should, in 
principle, provide a more reliable reference for the D/H determi- 
nation. 

The observing sequence foll owed the same schem e as that 
used for comet 103P/Hartley 2 (Hart ogh et all 1201 ll) . It con- 
sisted of ten 30-min long observations of the HDO 1 10— lot ro- 
tational line at 509.292 GHz, interleaved with 6-min simulta- 
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Table 1. HIFI observations of C/2009 PI (Garradd) on 6 October 201 1. Spectra characteristics and production rates. 
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(mK km s~') 
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HDO 


lio- 


-loi 


509.29242 


Single-point 


6.428 c 


22.5 ± 2.0 


-59 ± 33 


9.8(±0.9) 9.2(±0.9) 


7.7(±0.6) x 10 25 


HfO 


lio- 


■loi 


547.67644 


Single-point 


6.415 rf 


179 ±5 


-21 ± 13 


4.55(±0.12) 4.97(±0.12) 


3.94(±0.11) x 10 26 


H 2 


lio- 


-loi 


556.93600 


Single-point 


6.415^ 


9318 ± 12 


320 + 2 


23.8(+0.6) 23.4(±0.6) 


16.5(+0.4) x 10 28 


H 2 


lio- 


-loi 


556.93600 


Mapping 


6.291" 


8897 ± 98'' 


331 ± 20'' 


23.2(+1.6) 26.2(+4.0) 


20.8(±5.5) x 10 28 


H 2 


llO" 


-loi 


556.93600 


Mapping 


6.559 / 


8109 ± 82* 


293 ± 16* 


21.4(+1.8) 24.1(±3.7) 


19.2(±5.0) x 10 28 


H 2 Q 


2o2" 


-111 


987.92676 


Mapping 


6.597 g 


7140 ± 198* 


21 ± 19* 


18.2(+1.0) 20.6(±2.1) 


15.5(±2.8) x 10 28 



Notes. (a) Mean date. '*' Average of HRS and WBS line area retrievals in main-beam brightness-temperature scale T m % . The error bar corresponds 
to statistical noise, which is multiplied by 1.5 for single-point observations observed in FSW mode to account for uncertainties in the baseline 
removal. {c) Herschel observation identification number (Obsid) 1342230*, with # = 187, 189, 191, 193, 195, 197, 199, 201, 203, and 205. 

Obsid: # = 186, 188, 190, 192, 194, 196, 198, 200, 202, and 204. (c) Obsid: #185. (/) Obsid: #206. {s) Obsid: #209. w From the average of 
spectra within 10" of the peak brightness. (,) For single-point data, a 4" beam offset is considered (Appendix A). For mapping data, the production 
rate corresponds to the weighted mean of apparent production rates measured over the whole map, and the error bar is the dispersion around the 
mean. 
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Fig. 1. HIFI spectra of comet C/2009 PI (Garradd) observed on 6 
October 201 1 with the HRS. HDO (509 GHz), H 2 (557 GHz), 
and Hj 8 (548 GHz) lio-loi spectra are the average of single- 
point measurements. The spectrum of the H2O 2o2— 1 1 1 line at 
988 GHz is extracted from the map, by averaging data at offsets 
< 10" from the peak. The velocity scale is given with respect to 
the comet rest frame. Synthetic line profiles obtained with 30% 
extended production (see text) are shown by red dotted lines. 
Gas acceleration (velocity increasing from 0.48 km s _1 to 0.58 
km s -1 from r — 10 4 km to 10 5 km) is cons i dered to fit more 
closely the wings of the profiles (ICombi et al.Ll2004l) . 

neous measurements of the H2O and Hj 8 1 10— lot ortho lines 
at 556.936 GHz and 547.676 GHz, respectively. In addition to 
these single-point measurements, two on-the-fly maps of the 
H2O 1 10— lot transition (of 16-min duration) were acquired at 
the beginning and end of the sequence. Finally, a map of the 
H2O 2o2-l 11 para transition at 987.926 GHz was performed with 
an integration time of 25 min. The full sequence spanned 6.85- 
14.71 UT on 6 October, with t he maps serving to constrain th e 
H 2 excitation (lHartogh et al.Ll2010tlde Val-Borro et al.Ll2010h . 

The H2O and H^O 1 10— loi lines were observed in the upper 
and lower sidebands of the HIFI band la mixer, respectively. The 



HDO line was also observed with the band la mixer, whereas 
data on the H2O 2o2— 1 1 1 line at 988 GHz were acquired using 
the band 4a mixer. The single-point observations were carried 
out in the frequency-switching mode (FSW) with a frequency 
throw of 94.5 MHz. On-the-fly 557 GHz and 988 GHz maps 
were acquired using Nyquist sampling, and spatial coverages of 
4'x4' and 2'x2', respectively. The observing mode for the maps 
used a reference position at 10' from the comet in RA. Spectra 
were acquired with both the Wideband Spectrometer (WBS) and 
High Resolution Spectrometer (HRS). The spectral resolution 
of the WBS is 1.1 MHz. The HRS was used either in high- 
resolution (125 kHz) or nominal -resolution mode (250 kHz), en- 
abling us to sample the line shapes at a spectral resolution of 
70-150 m s . The telescope beam sizes at the frequencies of 
the three lines observed in band la are similar (half -power beam 
widths of 38'.'1, 38'.'7, and 4T/6 for the H z O, H' s O, and HDO 
lines, respectively), so that the three molecules were observed 
over the same (~ 55 000 km diameter) region of the coma. 

Figure Q] shows the HRS spectra of the single pointing mea- 
surements (HDO, H2O, and Hi 8 O 1 10— loi lines), as well as the 
H2O 2o2— 1 1 1 spectrum extracted from the central part of the 
map. Data reduction and calibration uncertainties are discussed 
in Appendix A. The HDO line is detected with a line-integrated 
signal-to-noise ratio of 17. Maps are shown in Fig [2] Line inten- 
sities and velocity offsets (Av) in the comet rest-frame are given 
in Table [1] The H2O 1 10— lot line at 557 GHz is optically thick 
and has an asymmetric profile owing to self-absorption in the 
foreground coma. Intensity variations with time of up to 7% are 
observed for this line, which is related to intrinsic comet variabil- 
ity or instrumental effects. Optically thin HDO and Hj 8 lines 
have approximately symmetric profiles. As the phase angle was 
32°, the small negative-velocity offset observed for these lines 
suggests a modest excess of outgassing toward the Sun. 

3. Analysis 

The analysis was carried out usin g one-dimension al excitation 
models of HDO, H 2 0, and H* 8 dBiver et al.Ll2007l) . Models in- 
clude collisions with H2O and electrons, which dominate the ex- 
citation in the inner coma, and solar infrared pumping of vibra- 
tional bands followed by spontaneous decay, which establishes 
fluorescence equilibrium in the outer coma. Radiation trapping 
strongly affects H2O excitation and is considered using the es- 
cape probability formalism. Consistent results were obtained 
using state-of-the-art radiation transfer methods (Hartogh et all 
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Fig. 2. On-the-fly maps of H 2 in C/2009 PI (Garradd) obtained 
with the WBS. Left: the 1 10 -lot 557 GHz line observed on 6.291 
October 2011 UT. Right: the 2 02 -lii line at 988 GHz observed 
on 6.597 October 201 1 UT. The contour spacing is 1 K km s _I 
in brightness temperature, corresponding to ~ 6cr. The Sun di- 
rection is indicated at lower right. The beam sizes are 38'.' 1 and 
21'.' 5 at 557 GHz and 988 GHz, respectively. 
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Fig. 3. Apparent H2O production rates as a function of beam 
offset deduced from the 557 GHz (left) and 988 GHz (right) 
H2O maps observed on 6.56 and 6.60 October UT, respectively. 
Calculations with a 47 K and variable 7\i n are shown with red 
and black symbols, respectively. The temperature profiles are 
plotted in the right panel with dashed lines. 



2010; Ben sch & BerginL l2004t IZakharov et alll2007h . Synthetic 
spectra are then computed using radiative transfer modelling. We 
assumed isotropic outgassing at a constant velocity. The model 
input parameters are: i) the gas expansion velocity, assumed to 
be 0.6 km s _1 , corresponding to the half widths of optically thin 
HDO, H 2 (988 GHz), and H^ 8 lines; ii) the gas temperature 
profile; iii) the x ne scaling fa ctor of the electron density pro- 
file, t aken to be equal to 0.2 dBiver et al.L 120071; iHartogh et ail 
2010; Ide Val-Borro et al.ll2010h . The hTq ortho-to-para ratio is 
assumed to be 3. 

Line intensities for spectra probing the inner coma are sen- 
sitive to the gas temperature profile, which controls both the 
population of the rotational levels and the optical depth of the 
lines. Hence, the evolution of the intensity of the H2O lio-loi 
and 2q2— 1 1 1 lines with beam offset, p (km), carries information 
about the temperature profile. Figure [3] presents the apparent 
H2O production rate 2app(H20) as a function of p deduced from 
the maps. With appropriate modelling, the Q app curve should 
be flat, if the nucleus is the dominant source of water vapour 
with a constant outgassing rate. Figure [3] presents the results 
for 7\i n = 47 K, a value consistent with multi-line observa- 
tions of methanol undertaken with the IRAM 30-m telescope 



in September and October 2011 wi th a 17" beam comparable 
to the HIFI beam dBiver et all 120121) . For both lines, g ap p(H 2 0) 
increases with increasing p for a constant temperature profile. 
Deviations from constan t Qann( H?Q) are enhanced when using 
x ne > 0.2 (lHartogh et al.L 120101) . The temperature law 7i aw that 
minimizes the deviations of 2 app (H20) from a constant value, 
and provides consistent (within 15%) water production retrievals 
from the two lines, has a minimum of 20 K at a distance r = 
4-20 x 10 3 km from the nucleus (thereby increasing the self- 
absorption and Qapp values) and then increases up to 150 K at r 
= 8 x 10 4 km (Fig. [3}. This temperature increase may be related 
to the increased efficiency of photolytic heating (ICombi et all 
l200l . We note that the line intensities are weakly sensitive to 
the temperature at r < 1000 km and r > a few 10 4 km. The ve- 
locity offsets Av of the on-nucleus synthetic 557(988) GHz line 
profiles are +292 (+50) and +259 (+39) m s -1 , for T law and T^n 
= 47 K, respectively. Hence, the model with a variable temper- 
ature also fits more closely the large positive Av of the central 
557 GHz spectra (+320 m s _1 , TableQ}, as it enhances the opti- 
cal thickness of the line. 

We examined whether the increase in 2app(H20) with p ob- 
tained with constant 7\j n could instead be attributed to water pro- 
duction from icy grains in the outer coma. For 7\m = 47 K, we 
are able to obtain a flat 2app(H20) profile for the 557 GHz line 
when assuming that 90% of water production is extended with a 
characteristic Haser scale-length of L ext = 30 000 km. However 
this model is unsatisfactory since: i) the apparent production rate 
given by the 988 GHz line now decreases with increasing p; 2) 
the predicted Av of the 557 GHz line is 30% lower than ob- 
served; 3) the inferred 16 0/ I8 ratio in water is then three times 
lower than the Earth value. On the other hand, subliming icy 
grains were possibly present in C/2009 PI (Garradd)'s coma 
(Pag anini et al.L 120121) . We therefore considered an alternative 
model with moderate (30%) water production from long-lived 
icy grains (model outputs are insensitive to an unresolved source 
of water). This model explains the water 557 and 987 GHz maps 
for L ext = 50 000 km and a temperature law similar to T\ avj (min- 
imum of 20 K at r — 4-16 x 10 3 km). The retrieved production 
rates and production rate ratios are identical (within 3%) to those 
found for nuclear production with T\ m , This model accounts sat- 
isfactorily for the observed line profiles (Fig. [TJ. 

Table Q] presents the production rates calculated with T\ m , 
along with those for constant T^ B = 25 K and 47 K. To compute 
the error bars in the production rate and isotopic ratios (derived 
from the simultaneous single-point data), we take into account 
a 5% relative calibration uncertainty (Appendix A). Using Ji aw , 
the HDO/H^O production rate ratio is 0.215 ± 0.023. This is 
significantly lar ger than the value .161 ± 0.017 measured for 
103P/Hartley 2 dHartogh et alll201 ll) . For 7\i„ = 25 K and 47 K, 
one finds HDO/H^O = 0.185 and 0.195 (+0.020), respectively, 
hence the retrievals are only slightly temperature-dependent. 

The H^O/H^O production rate ratios derived from the band 
la observations using T^ n = 47 K is 419 + 26, while for T\- d „ and 
Tkin = 25 K one finds 523 + 32 and 470 ± 29, respectively. These 
latter valu es are in good agre ement with previous measurements 
in comets (Jehi n et all 120091 and references therein), as well as 
with the 16 0/ 18 = 498.7 VSMOW value, giving further con- 
fidence that our radiation-transfer model properly accounts for 
the opacity of the Hi 6 O lines, provided a low gas temperature is 
adopted. 

We adopt the model with Ti aw for the determination of the 
D/H ratio, as it best explains the H2O maps and line profiles. 
The D/H measure ment for 103P/Hartle y 2 was based on the 
HDO/H' 8 ratio dHartogh et all 1201 ll) . assuming 16 0/ 18 = 
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500 + 50. Using the same method, we derive D/H = (2.15 ± 
0.32) x 10" 4 for C/2009 PI (Garradd). Using instead directly 
the HDO/H 2 16 production-rate ratio results in D/H = (2.06 ± 
0.22) x 10~ 4 . This value will be adopted for the discussion. We 
note that consistent results are derived using constant (25—47 K) 
gas temperatures and 30% extended production (central values 
ranging from 1.96 to 2.33 x 10~ 4 ) . 

The water production rate derived using T\. dvJ is in the high 
range of values measured by other techniques in October 2011 
(see discussion in Appendix B). Using a lower value would im- 
ply that the water of C/2009 PI (Garradd) is highly enriched in 
18 relative t o both the Sun and all rocky bodies of the inner 
solar system (McKeegan et ail 1201 lh . On the other hand, only 
an extreme ls O enrichment ( 16 0/ 18 ratio of ~ 330, i.e., <5 ls O 
= +500 %o in geochemical notation) would reconcile the D/H 
ratio in comet C/2009 PI (Garradd) with the canonical Oort- 
cloud value of 3 x 10~ 4 . Such high 18 enri chments have not 
been found so far in any solar system body dMcKeeean et all 
), except for CO in the Titan's atmosphere dCourtin et all 
Our preferred interpretation is thus that the 16 0/ 18 ratio 



2011 



2011 



in comet Garradd is consistent with the VSMOW value. 

4. Discussion 

The discovery of a D/H value equal to that of the Earth's oceans 
in the Jupiter-family comet 103P/Hartley 2 showed that the 
reservoir of Earth-like water in the solar system is substantially 
larger than previously th ought, including now both carbonaceous 
meteorites and comets (Hartogh et all l20llh . It also revealed 
that isotopic diversity is present in the comet population, with 
members of the Oort cloud having a deuterium enrichment of 
up to a fac tor of two with respe ct to VSMOW (Fig. As dis- 
cussed by Hartogh et all d201 lh . the suggested dichotomy be- 
tween Oort-cloud and Jupiter-family comets is difficult to ex- 
plain in the context of current models predicting deuterium en- 
richments in cometary ices, and, if real, would imply to revisit 
the source regions of OCCs and JFCs. Actually, the only dynam- 
ical theory that could explain in principle an isotopic dichotomy 
in D/H is the one arguing that a substantial fraction of Oort-cloud 
comets were captured from other stars when the Sun was in its 
birth cluster dLevison et alll2010l) . 

The paradigm for a single, archetypal D/H ratio for all 
Oort-cloud comets is no longer tenable. The value of (2.06 
+ 0.22) x 10~ 4 measured in comet C/2009 PI (Garradd) is 
smaller than the mean of previous determination s in O ort-cloud 
comets ((2.96 + 0.25) x 10~ 4 , lHartogh et all l201ll) . and is 
consistent with the upper limit of 2.5 x 10~ 4 m easured in 
com et 153P/Ikey a -Zhan g (Fig. gj iBiver et all l2006h . We note 
that iBrown et alj d2012f) reexamined mass-spec trometer mea- 
surements for lP/Halley dEberhardt et all 1 19951) . reevaluating 
this value to be 2.1 x 10~ 4 . Altogether, the available data suggest 
that the deuterium enrichment in the water of Oort-cloud comets 
is not as high as previously thought, at least for a fraction of 
the population. Nevertheless, the D/H ratio measured in comet 
C/2009 PI (Garradd) is significantly higher (by 2.3 rx) than the 
VSMOW value measured in 103P/Hartley 2. Interestingly, the 
range of D/ H ratios measured in St ardust samples from comet 
8 IP/Wild 2 dMcKeegan et all 120061) is the same as measured in 
the bulk water of various comets. 

Dynamical modelling suggests that the distribution of plan- 
etesimals underwent large-scale mixin g during the stages of 
planetary migration (Walsh et all 1201 il) . Therefore, variations 
in the cometary D/H may be expected within each population, 
if ancestor reservoirs were isotopically different. Isotopic diver- 
sity may be linked to the large compositional diversity observed 
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Fig. 4. D/H ratio in the w ater of come t s (inc luding the re- 
vised value for Halley from Brown et Id] d2012l) ) compared to 
values in carbonaceous meteorites (CI), the Earth's oceans, 
and Enceladus. Displayed data for planets, interstellar medium 
(ISM), and th e protosolar nebula refer to the value in H 2 . 
Adapted from lHartogh et all d201ll) . 

withi n both OCC and JFC populations (Bockele e-Morvarl 
1201 ll) . Finally, experimental studies of ice sublimation suggest 
that the D/H measured in the evaporated vapour might be en- 
hanced o r depleted with resp ect to the bulk D/H in the cometary 
nucleus dBrown et all 120121) . This demonstrates the need to in- 
crease the sample of comets with accurate measurements of the 
D/H ratio, as well as perform further modelling. 
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Appendix A: Data reduction and HIFI calibration 

The comet was tracked using an up-to-date ephemeris provided 
by the JPL Horizons system. The Herschel rms pointing accu- 
racy is approximately 1". 

The data were reduced to level 2 products using the Herschel 
Interactive Processing Environment (HIPE 7.3). All lines were 
observed in the two orthogonal H and V polarizations. The two 
orthogonal polarizations were averaged. Note that the two polar- 
izations are observed with different mixers, and their respective 
apertures are imperfectly co-aligned. The beam offset for the H 
and V average spectra is ~ 4" with respect to the pointed posi- 
tion. 

The line intensities integrated over velocity were computed 
on the main-beam brightness-temperature scale using beam effi- 
ciencies of 0.75 and 0.73 for bands la and 4a, respectively, and a 
forward efficiency of 0. 96. Based on the calibration error budget 
dRoelfsema et al.L [2012b . a conservative value for the uncertainty 
in the absolute intensity calibration is 10% for both bands. Most 
sources of errors are eliminated when comparing band la data, 
and the relative uncertainty is at most 5% in this case (sideband 
ratio, hot-load coupling, and temperature). Finally, the calibra- 
tion uncertainty for the ratio of band la to band 4a lines is 10%. 
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Appendix B: Water production rate: comparison 
with other measurements 

Water production rates measured for C/2009 PI (Garradd) from 
September to October 201 1 are shown in Fig. [5] Reported val- 
ues include retrievals from OH 18-cm observations (Colom 
et al. 2011, and in preparation), OH narrowband photometry 
(Schleicher et al., personal communication), and near-IR ob- 
servations of water (DiSanti et al. 2012, Paganini et al. 2012; 
Villanueva et al. 2012). Herschel/HIFI and OH 18-cm observa- 
tions, which were acquired close in date, provide consistent val- 
ues. The data suggest that the activity of comet Garradd under- 
went significant variations, and reached a maximum at the time 
the HIFI observations were performed. The 987 GHz H2O line 
observed with HIFI was also observed on 16 October 201 1 with 
the Spectral and Photometric Imaging Receiver (SPIRE) aboard 
Herschel; the line intensity is ~ 40% weaker than on 6 October 
2011, when the HIFI measurements were conducted (Swinyard 
et al., in preparation), implying a water production rate consis- 
tent with the value measured from OH narrowband photometry 
on 18 and 20 October 2011 (65 days before perihelion, Fig. [5). 
The low values derived from the near-IR measurements, com- 
pared to other measurements, possibly reflect sublimation from 
short-lived icy grains in the inner coma since the field-of-view 
for these observations was much smaller (by a factor of ten or 
more) than for the other techniques. 



Fig. 5. Water production rates measured pre-perihelion in comet 
C/2009 PI (Garradd). The time is with respect to the perihelion 
(23 December 2011). The plain and empty circles correspond 
to the 6-October single-point and mapping HIFI observations, 
respectively, analysed with the model with 7\ aw . 



